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ABSTRACT. This article provides sufficient conditions for positive maps on the Schatten classes 
Jv> 1 < p < oo of bounded operators on a separable Hilbert space such that a corresponding 
Perron-Frobenius theorem holds. With applications in quantum information theory in mind suf- 
ficient conditions are given for a trace preserving, positive map on J\, the space of trace class 
operators, to have a unique, strictly positive density matrix which is left invariant under the map. 

f^*) , Conversely to any given strictly positive density matrix there are trace preserving, positive maps 

f*"^ ' for which the density matrix is the unique Perron-Frobenius vector. 

o 

(N 

Dedicated to S. Doplicher and J.E. Roberts on the occasion of their 60th birthday 

1. Introduction 

In the theory of quantum information the transmission through noisy channels plays an im- 
portant role. Usually it is described by what physicists either call a quantum operation (see e.g. 



X 



[25]) or a stochastic map ([fil], see also [|i"9|]) or a super-operator (see e.g. [28]) and what mathe- 



maticians call a completely positive (trace preserving) map. In mathematics positive maps were 



£ — . first studied by Kadison in the context of C*-algebras [ ]18[ ] and completely positive maps by 

O Stinespring [ |33|]. In quantum physics it was first introduced by Haag and Kastler, who called it 

an operation []15|]. They were then studied in more detail in [16]. [21, 



o 

Oh 



contains an extensive 



discussion, how this concept naturally arises in quantum physics. 

Because they preserve positivity such maps naturally fit into a context where it makes sense 
to ask about a formulation of a corresponding Perron-Frobenius theory. 
| As the name indicates, the earliest results relating positivity of a finite dimensional linear map 

to the non-degeneracy of an eigenvalue and the positivity of the resulting (Perron-Frobenius) 
eigenvector is due to Perron and Frobenius [ p7[ , |l S]. The first extension to the infinite dimen- 
sional was given by Jentzsch JT7|]. Since then there has been an extensive development (see e.g. 



, 23, 30]) 



The interest in quantum physics and quantum field theory originates from the observation, 
that very often the ground state of a quantum system with Hamilton operator H > is non- 
degenerate and nowhere vanishing. Indeed, under suitable circumstances such a ground state 
can be viewed as the Perron-Frobenius vector for exp —tH, t > 0, an observation first made by 
Glimm and Jaffe [||] (see also ^ ||] for an overview and with references to other articles). 

The Perron-Frobenius theory and related topics for positive linear maps on linear spaces of 
operators was first studied in [f7], 14, |0, where in the first article the space of operators 



was the C*-algebra of all linear operators on a finite dimensional Hilbert space. In the second 
and third articles the discussion was extended to the infinite dimensional case and in the fourth 
article the analysis was carried out in the context of a von Neumann algebra. 

The present article can be viewed as an extension of the discussion in [^] to the infinite 
dimensional case. Indeed, we will consider positive maps on the Schatten classes J p , 1 < p < 
oo and provide a Perron-Frobenius theory on such spaces. As in the usual context the notions 
of ergodicity and positivity improving will play an essential role. Since density matrices are 
elements of J\ our main focus will be on the case p = 1. In particular we will show that 
any compact, ergodic, completely positive and trace preserving map in J\ has a unique density 
matrix invariant under this map (i.e. the eigenvalue equals 1) and which therefore can be viewed 
as the Perron-Frobenius vector for this map. 
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The article is organized as follows. In Section 2 we will briefly review the concepts needed 
for our discussion. In Section 3 we will prove a Perron-Frobenius theorem (non-degeneracy of 
a certain eigenvalue), provided such an eigenvalue exists, a typical condition needed in infinite 
dimensional contexts. Section 4 will provide sufficient conditions for such an eigenvalue to 
exist. In Section 5 we will give examples. Also we will prove a converse in the sense that to 
any density matrix > there are completely positive maps for which this density matrix is the 
Perron-Frobenius eigenvector. 



2. Positive and Completely Positive Maps 

We start by recalling and establishing some definitions and facts related to trace ideals and to 
the concept of positive and completely positive maps ( for more details on trace ideals see e.g. 



03U |32Q and on positive and completely positive maps e.g. |f33j, |34j, |5|, g, |35|, |26|]). Let H be any 
complex, separable Hilbert space with scalar product denoted by ( , ). By B(TL) we denote the 
C*-algebra of all bounded operators on TL equipped with the norm 1 1 1 1. By I € B(TL) we denote 
the identity map on TL. 

For any A G B(TL) we set \A\ = (A* A) 1 / 2 G B(TL) where * denotes the adjoint. Also we 
let Tr be the trace operation on TL. For 1 < p < oo let J p = J P (TL) be the Schatten class 
consisting of all elements A G B(TL) such that \A\ P is trace class. J p is equipped with the 
norm ||-A|| P = (TV(| J 4| P )) 1/ ' P making each J p a complete Banach space. One has A* G J p 
if A G J p , i.e. the space J p is self-adjoint in B(Ti), such that ||^4*|| p = Also we set 

Joo = B(TL) again with the same norm, i.e. || = || ||. We will frequently make use of the 
trivial identity \\A\\ P = \ \ \ A\ \ \ p valid for all A G J p and all 1 < p < oo. Via the trace J q is 
the dual of J p with l/p + 1/q = 1 for all 1 < p < oo, i. e. each continuous linear functional 
on J p is given by an A G J q in the the form Tr(A*B), B G J p . In particular the spaces J p 
are reflexive for all 1 < p < oo. We prefer to use the adjoint in describing linear functionals 
since for p = 2, in which case J2 is the space of Hilbert-Schmidt operators on TL, the scalar 
product making it a Hilbert space is just given by (^4, B)2 = Tv(A*B). Also one has the Holder 
inequality | Tr(A*B)\ < ||A|| 9 ||5||p. Let C = C(TL) = Coo be the closed positive cone of all 
elements A > in B(Ti). As usual we write A > and A > if the relations (<p, Aip) > 
and (ip, Aip) > respectively hold for all 7^ (p G TL. Correspondingly we write A > B 
(or B < A) and A > B (or B < A) if A - B > and A - B > respectively. Obviously 
A>B>0orA>B>0 implies A > 0. As is common, A is said to be positive definite 
if A > 0. Set C p = C P (H) = C n J p , a closed set in J p . Also A > in J p if and only if 
Tr(BA) > for all < B G J g . Correspondingly A > in J p if and only if Tr(BA) > for 
all < B ^ in J q (l/p + 1/q = 1, 1 < p < 00). In this sense the cones C p and C q are dual 
to each other. The closed set Ci t i = {A G C\ \ Tr A = \\A\\\ = 1} in J\ is the set of all density 
matrices in TL. 

By definition a positive map <p in J p is a linear map from J\ into itself, which leaves C p 
invariant. 4> is called n-positive if the induced map cj) n = (g) I n in Jp (g) B(TL n ) also leaves the 
corresponding cone C P (TL <8> TL n ) of non-negative elements invariant. Here TL n is any Hilbert 
space of dimension 1 < n < 00. Obviously if 4> is n-positive, then it is also n'-positive for any 
1 < n' < 11 as is A0, A > 0. If cj> and <p f are both n-positive on J p then so is their composite 
(f> o (f>' and their sum <p + <fi'. So the n-positive maps in J p form a cone in the linear space of all 
linear maps in J p . 

If <f) is n-positive for all n, then is called completely positive. Thus the map (j) a : A — > 
aAa* for any a G B(TL) is completely positive on J p for all 1 < p < 00 as is any finite linear 
combination cf) a = J2i <Pcti with a-i G B(TL). If a has an inverse a -1 in B(TL) then cf)^ 1 = <p a -i. 
Also cf) a o (j) a i = (j) aa / and 4>\ a = |A| 2 (/> a for a, a' G B(TL) and A G C. As it turns out at least 
for p = 1 the induced linear map <j) <£> Aw from J p= \ <&B(TL') then also leaves the corresponding 
cone invariant for any separable Hilbert space TL 1 (see e.g. [f22|]). Since density matrices are in 
Ji, the case p = 1 is of most interest in quantum physics. Then there are good physical reasons 
to consider completely positive maps rather than only positive maps (see e.g, [f22|]). 

Although the following lemma is well known in similar contexts, we still will provide the 
short proof. 
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Lemma 2.1. The following relation holds if the map cp in Jp (1 < p < oo) is positive. 

(1) 0(A)* = 0(A*), A G J p . 

Proof. We first consider the case when A is self-adjoint, i.e. A* = A. Write A = A + — A_ 
with A± E C p , A+A_ = A_A+ = and \A\ = A + + A_. Here ±A± are the positive and 
negative parts of A, obtainable from the spectral representation of A or more explicitly as A± = 
l/2(\A\ ±A). Obviously \\A±\\ P < \ \A\\ p . Since <j>{A±) > 0, 0(A) = cp(A + ) - cp{AJ) is self- 
adjoint. For arbitrary A G B{H) write A = ^RA+iQA G B{H) with KA = 1/2(A+A*), OA = 
l/2i(A - A"), such that KA, OA E B(H) are self-adjoint with KA, OA G J p whenever 
A G Jp. More precisely we have the a priori bound ||3?A|| p < ||A|| p , ||9A|| p < ||A|| P . This 
gives the decomposition 

A = (KA) + - (5ftA)_ + *((9fA) + - (3A)_) 
and hence by the linearity of 

(2) 0(A) = 0((3?A) + ) - 0((KA)_) + *0((3A) + ) - i0((3A)_) 

with ||^(^)|| p < ||0((^) + )|| / ,+||0((^)_)||^+||0((^) + )|| p +||0((^)_)|| ?3 . Inparticular 
(Q) shows that indeed (g) holds. □ 

We also note the following. The relation -\A\ < A < \A\ for A = A* € J p implies -0(|A|) < 
0(A) < 0(|A|), whenever the map cp in Jp is positive. This in turn gives 

(3) |0(A)|<0(|A|) 

valid for any i = i* £ J p and any positive map in JL. We do not know whether (or 



when) this relation continues to hold when the condition A = A* is dropped (see, however, < J23[ ) 
below for a weaker result, which will suffice for our purposes). Observe by comparison that for 
any n x n matrix 5 with non-negative entries, the classical context for the Perron-Frobenius 
theorem, one has \Sz_\ < S\z} for z G C n . Here \z_\ is the vector whose components are the 
absolute values of the corresponding components of z. Also for real vectors x = {x{\ and 
y = {Hi} by definition x < y if and only if xi < yi holds for all i. 

Furthermore one does not necessarily have cp(A)± = 0(A±) for any A = A* G J v . 
However, since — A_ < A < A + implies — 0(A_) < 0(A) < 0(A + ), the inequalities 
0(A)± < H A ±) hold - 

Lemma 2.2. Any positive map cp in Jp (1 < p < oo) is continuous, i.e. it satisfies \ \cp\\ p < oo. 
Here we denote by | \cp\ \ p the norm of any continuous linear map cp in J p , such that 

(4) \\cp\\ p = sup | || p . 

I|A|| P <1 

In particular 

(5) ||0||i = sup TY(0(A)), 

AeCi,i 

if cp is completely positive on J\ (see e.g. [|6|, |22j] ). 



Proof. We adapt a standard proof (see e.g. []26|], p. 19) used in the context of positive maps on 
non-unital C7*-algebras. We first claim that it suffices to prove boundedness on C p . Indeed, this 
follows from the decomposition @ and the related bounds. Assume that cp is not bounded on 
C p . Then there are A n G C p with ||A„|| p < 1 and ||0(A„)|| P > n 3 . Let A = J2 n l/"- 2 A„ G C p , 
such that < l/n 2 A n < A holds for all n. 
We need the fact that 

(6) \\B\\ P <\\B'\\ P 

holds for any < B < B' G C p . The case p = oo is trivial. When 1 < p < oo take ip n to be a 
complete orthonormal basis in TL diagonalizing B. Then we have 

\\B\\l = ^2(<p n ,Bcp n f < ^2(cp n ,B'cp n )P 

n n 
< ^(<p n ,B ,P <p n ) = \\B% 
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proving (j^). Here we have used the estimate 

(7) \{ip,A<p)\P<{<p,\A\Pip), 



valid for any A = A* S B(H), any normalized ip € 7i and any 1 < p (see e.g. [|32|], p.21). 
Since is positive (§) gives n < l/n 2 \ \(f)(A n )\ \ p < \ \(f)(A)\ \ p , which is a contradiction. □ 

For any continuous linear map 4> in J P , 1 < V < oo let <fi* be the "adjoint" continuous linear 
map from J q , l/q + l/p = 1 into itself given by the relation Tr(((j)*(A))*B) = Tr(A*4>(B)) 
for all A € J q and B G J"p. In particular 0** = holds whenever 1 < p < oo. Also 0* = a * 
with a* = (ai, ...a*, ...) by the cyclicity of the trace. 

By definition a continuous linear map <ft in J\ is called trace preserving if Tr((j)(A)) = Tr(^4) 
holds for all A G J\ and this in turn is equivalent to the relation <fi* (I) = I. Obviously trace 
preserving maps leave the set C\ } \ of density matrices invariant. If (j) in J\ is completely positive 
and trace invariant, then by (0) \ \4>\\i = 1- 



All completely positive maps on J\ have the Kraus representation [22], a consequence of a 
theorem of Stinespring [^] (see also [Q] for a proof in the finite dimensional case): Given such 
a cj) there is an at most denumerable set of elements a = {aj}i g N in B(7i) satisfying 

(8) £Xoi< ||0| |il 



for any finite subset K C N such that = 0^, again with 

(9) a (A)=^aj^, 



which now may be an infinite sum. This representation is not unique. Conversely each such 4> a 

is completely positive. More precisely (j) a is defined as follows. Let a N = (a\, , «at). Then 

the cj) a form a Cauchy sequence with respect to the norm [| ||i and (p a is defined as the limit 



(see e.g. [22]). If J2i ct*a>i = I, such that ||0||i = 1, then 4> a is trace preserving in J\. In the 
finite dimensional case (dim'H < oo) the representation may always be chosen such that the 
index i runs through a finite set of the order at most (dim TL) 2 . 

An important and interesting feature of completely positive maps in the context of quantum 
physics is that they not necessarily map pure states, i.e. one-dimensional orthogonal projections, 
into pure states. 

Again we have the following representation for the adjoint of 4> a 

(10) ^(A)=J2a*A ai , 

i.e. 0* = 4>a*. So 4>a is trace preserving if and only if 0^,(1) = I an d then | |0qJ|i = 1. Observe 
that to any completely positive (ft with \ \4>\\i < 1 or equivalently 0*(I) < I ( this condition 
is again natural in the context of quantum physics, see e.g. [^2[]) we may in a natural way 
associate a trace preserving, completely positive map <p given as 4>{A) = 4>{A) + aAa where 
a = (I-0*(I))V2 > o. 

An example to which we shall return below is when the a, 's are orthogonal projection oper- 
ators which are pairwise orthogonal, i.e. satisfy 

(11) ajQj = SijCti = 5ija*. 

Then <jf a (J) = I if and only if the aj provide a decomposition of unity, i.e. if J2ien a i = ^ 
holds. As already remarked (p a is then also trace preserving. More generally if the a-i are all 
selfadjoint with a i = I> then <ft a is trace preserving. 

We remark that there are some situations (see e.g. j pO| ] and which actually was the motivation 
for the present discussion) where one starts with completely positive maps of the following form. 
Let (0, n) be a measure space, i.e. /x is countably additive. Suppose we are given a measurable 
map uj — > a(uj) from £1 into B{TL) such that 

(12) / \\a{uj)\\ 2 dii(uo) < oo. 
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We then set 

(13) <j)(A) = J a(Lu)Aa(Lu)*dfi(cj), 
such that its adjoint takes the form 

(14) <f>*(A) = J a(io)*Aa(Lo)dfj,(Lo), 
The condition corresponding to (|J) is given as 

(15) <j>*(I) = [ a(uj)*a{u))dn{u)) < cl 



for some < c < oo. To cast <p in the form (g) such that ([TJ|) gives (||), let Xi be any orthonormal 
basis in the separable Hilbert space L 2 (Q, dfi) and set 



Oii = J Xi(u)a(uj)dii(u). 

As an example one may choose SI to be a compact group, /j, its Haar measure and u — > a(ui) a 
unitary representation. Then the image of <fi consists of all trace class operators which commute 
with all a(u>). Also <fi is trace preserving and an idempotent, i.e. <p 2 = <p. Moreover <p is a 
conditional expectation when dimTC < oo, i.e. <fi(A(j)(B)) = <p(A)<p(B) = 4>{4>{A)B) holds 
for each A,B. 

3. The Perron-Frobenius Theorem 

Before we turn to a discussion of the Perron-Frobenius theorem we make some general re- 
marks on (J p (4>) for the map 4> in J p , in particular when (j) is positive. Here a p (4>) C C is the 
spectrum of (j), i.e. the set of all A for which (A — <p) does not have a bounded inverse in J p . Let 
r p {4>) be the spectral radius of any bounded linear map <p in J v , i.e. 

(16) r p (<f>) = lim U n \\ l J n < 



By a celebrated general result of Gelfand Q11Q, one has sup |c p (</>)! = r p{4 ) )- If the map (f> in 
J\ is trace preserving and completely positive then r\{&) = ||(/>||i = 1. Indeed, 4> n is then also 
completely positive and trace preserving such that by (|J) | |0 n | |i = 1 holds for all n. 

Lemma 2.1 allows us to draw the following conclusions on a p {4>) for positive maps in J p . 
First we observe that in general a p (<p) is not contained in the real axis. In fact, in the finite 
dimensional case a p ((j)) C R if and only if = <p and in the infinite dimensional case the 
same statement is valid when p = 2. Let p p {4>) be the resolvent set of the map 4> in J p , i.e. the 
complement in C of a p {4>). If A € p p {4>), then to each A G J p there is a unique A' £ J p , such 
that (X-(j))(A') = A. Taking the adjoint and using ([]]) gives (X-(p)(A'*) = A*. Since A can be 
chosen arbitrary also A* can be chosen arbitrary giving a unique A* , this shows that both sets 
p p {4>) and cr p {4>) lie symmetric with respect to the real axis. Let A G cr(<^>) be a real eigenvalue 
with eigenvector A: <p(A) = A A. So by ([]]) A* is also an eigenvector with the same eigenvalue. 
With the decomposition A = A\ + iA2 we see that both self-adjoint operators Ai and A2 are 
eigenvectors, so any eigenspace for a real eigenvalue is spanned by self-adjoint elements. If in 
addition the map cj> in J\ is trace preserving and A / 1 an eigenvalue and A a corresponding 
eigenvector, then from Tr(^4) = Tr(^(^4)) = A Tr(A) we deduce Tr(A) = 0. 

In the finite dimensional case dim'H < 00 all spaces J p are equal and of course a p {4>) is 
independent of p. Since all norms 1 1 | \ p are equivalent, by its definition r p {4>) is also independent 
of p, as it should be. 

To formulate the Perron-Frobenius theorem in the present context, we make the following 
definition, which is the just an adaption of the usual definition. Since any positive map (f> in J p 
is bounded exp t<j) is a well defined positive map in J p for all t > 0. Also its inverse exp — t<p is 
a well defined bounded map in J p , not necessarily positive . 

Definition 3.1. A positive map <fi in J p (1 < p < 00) is positivity improving if<j>(A) > Ofor 
any A > 0, A 7^ 0. (p is ergodic if for any A > 0, A ^ there is tA > with (exp tA<f>){A) > 
0. 
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If is ergodic then (exp t(f))(A) > (exp tA<p){A) > for all t > tA- Obviously is ergodic 
if it is positivity improving. 

A simple necessary criterion for ergodicity is given by 

Lemma 3.1. If<f> is ergodic then Kercj)\c p = 0. 

Proof. Assume that there is < A / with <fi(A) = 0, such that (exptcj))(A) = A for all 
t. If A is not positive definite, then we have a contradiction, so let A > 0. We claim there is 
< A' / which is not positive definite with A' < A. In fact, we may take A' to be a one- 
dimensional orthogonal projection onto an eigenvector of A times the corresponding eigenvalue, 
which we choose not to be zero. This gives < 4>{A') < <f){A) = and we are back to the first 
situation. □ 

Some remarks concerning the finite dimensional case are in order. 

First we will show that when dim TL < oo the present definition of ergodicity of is equiva- 
lent to the condition (1 + 0) dim ' H ~ 1 (,4) > for any < A / 0, a criterion used in [^]. Obvi- 
ously this last condition implies ergodicity since expia > l/(n!) 2 min(l, + a) n for all n 
and all t, a > 0. As for the converse assume there is < A ^ such that (1 + <j)) dlvaH ~ 1 {A) 
is not positive definite. So there is an orthogonal projection P ^ 0, 1 and A > such that 
0(-P) < AP (see [Q]). This gives (expt0)(P) < (expiA)P for all t > 0, so is not ergodic. 

We have the following necessary and simple criterion for to be ergodic. 

Lemma 3.2. If dimTL < oo then (f)(1) > Ofor ergodic 0. 

The converse is not true as may be seen by taking to be given by @, where the Oj's are 
taken to be a nontrivial decomposition of unity. Then 0(1) = I whereas 0(a«) = ai for all i, 
such that 4> is not ergodic. 

Proof. Assume that 4>(T) > holds such that there is ip E TL with (tp,<j)(T)(p) = 0. Then 
(ip, (j)(A)(p) = holds for all selfadjoint A due to the bound -| \A\ |I < A < \ \A\ \I and therefore 
(<p, (f) n (A)ip) = for all n > 1 and all A = A*. This in turn implies (p, (exp t(j))((j)(A))<p) = 
for all i > and all A = A*. Since ^ by definition, there is A > with < 0(A) / (see 
the decomposition (0)), so <p is not ergodic. □ 



Returning to the general case it is clear that a given by (g) is not ergodic if the aj 's satisfy 



condition (|llj). Also in the context p = 1 and with a continuous, one parameter family of 
unitary operators U(t) defining a quantum mechanics on TL in terms of its infinitesimal generator 
H, then 4>u{t) (which describes the time evolution in the Heisenberg picture) for fixed t is 
not ergodic. If the completely positive map in J v is given by (13) in terms of a unitary 



representation of a compact group, which contains at least one irreducible representation with 
finite multiplicity, then by Schur's Lemma is easily seen to be ergodic ( and even positivity 
improving) if and only if the representation is irreducible. Then also TL is of course finite 
dimensional and the map is a conditional expectation. The following two lemmas are almost 
obvious. 

Lemma 3.3. A completely positive map in J\ represented as <p a is positivity improving if and 
only for any ^ <f> € TL the closed linear hull of the set of vectors {a*(p} is all ofTL. 

Proof. It suffices to consider 4>{P), where P is any one-dimensional orthogonal projection 
(compare the proof of Lemma 3.1). Let 7^ <po G PTL be a unit vector, so we have to consider 

(17) (<p,<KP)<p) = El^ a ^o)! 2 

iGN 

and the claim follows. □ 

This proof shows that in the infinite dimensional case no a with only finitely many non- 
vanishing ttj's can be positivity improving. In fact, for given P choose 7^ ip to be orthogonal 



to all aji^o- Then (17) vanishes and therefore 4>{P) is not positive definite. 



Let A(a) be the algebra (not *-algebra) generated by the a*'s. 



PERRON-FROBENIUS THEORY FOR POSITIVE MAPS ON TRACE IDEALS 



7 



Lemma 3.4. A completely positive map <f> in J\ represented as (p a is ergodic if and only if every 
non-zero vector ip in Ji is cyclic for -4(a), i.e. the strong closure of the linear space A(a)(p is 
all ofTL. 

These two results apply to all Kraus representations 4> a of 4>. 

Proof. We use the same notation as in the proof of the previous lemma. In addition write 

a, = a il ai2-.a in for i = (i\,i2, -i n ) and set |i| = n. Then 

oo , n 

(tp, (expt<f>)(P)(p) = -j K^Oi^o)! 2 

n=0 ' i:|i|=ra 

and again the claim follows. □ 

The next theorem is part of the Perron-Frobenius theorem in the present context. It is part 
because it assumes that r p (<p) is an eigenvalue, a standard assumption one has to make in the 
infinite dimensional context if no other information is available. 

Theorem 3.1. Let the positive map (ft in J p (1 < p < oo) be ergodic and assume that r p ((j)) = 
\\4>\\ p . Ifr p ((f>) is an eigenvalue of (j), then this eigenvalue is simple and the eigenvector A may 
be chosen to be positive definite. 



Remark 3.1. According to usual terminology such a vector is called a Perron- Frobenius vector. 
When dm\7i < oo then r p ((fi) is independent of p and is always an eigenvalue. This case has 
been covered by Evans and Hoegh-Krohn [f7|] with an application to quantum Monte-Carlo 
processes and by Groh [ 14]. At present we do not know whether in the infinite dimensional case 
the condition r p ((f)) = \\4>\\ p can be dropped. 

For the proof of the theorem we need the following simple 

Lemma 3.5. For any B > and ^ B' > in J p , 1 < p < oo the following strict inequality 

(18) \\B-B'\\ P <\\B + B'\\ P 
holds. 

The condition B > cannot be weakened to B > as the choice B = F / and B' = 
Q ^ for two finite dimensional orthogonal projections with PQ = shows. 

Proof. The estimate ( |l8| ) for the case p = oo is trivial, so let p < oo. Since B — B' € J p 
is self-adjoint, there is a complete set orthonormal eigenvectors tp n with real eigenvalues A n of 
B - B'. With A n = (ip n , (B - B')<p n ) we obtain 

(19) \\B - B% = |An| P = £ \(<Pn, Bifn) - (<p n , B' ^ . 

n n 

Observe that 

\{cp n ,Bip n ) - {(p n ,B'ip n )\ p < m&yL(((p n ,Bp„) p ,(p n ,B'ip n ) p ) 

< (((p n , Bip n ) + {(f n , B'ip n )) p 

= {ip n ,(B + B')<p n )P 

(20) < (tp n ,(B + B'f(p n ) 

holds, since B > and B' > 0. In the last estimate in ( j20| ) we used the estimate (0). 

Furthermore (ip n ,Bip n ) > holds for all n. Also (ip n ,B'ip n ) > for at least one n. 
Otherwise we would also have {tp n , B'ip n i) = for all n, n' by Schwarz inequality and hence 
B' = 0, contradicting the assumption. Thus the first inequality in (j20|) is actually strict for at 



least one n. Inserting (2C) into (19) and using this last observation proves (18). □ 
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Proof of the Theorem: Let A ^ be any eigenvector for <p with eigenvalue r p (<p). By an 
observation made above, we may assume that A is self-adjoint. Observe the following estimate, 
valid for any t > 0. 

exptr p (4>)\\A\\ p = \\exptr p (^)A\\ p = \\(expttp)(A)\\ p 

= |||(exp^)(A)||| p < ||(expt0)(|A|)|| p 

(21) < expt||^||p|||A|||p = expt||^||p||A|| p . 

Here we have used the estimate (§) in combination with (|3) for the map exp tcft. 

Since by assumption r p (<p) = \ \<f>\ \ p , all inequalities in (2l]) actually have to be equalities. In 
particular for all t > 

1 1 (exp t(f))(A)\\ p = || (exp t<f>) ( A + ) - (exp t<j>) 1 1„ 

(22) = ||(expt0)(|>l|)||p = ||(exp^)(A + ) + (exp^)(A_)||p. 

If A + = or A_ = there is nothing to prove so assume that A± / 0. By the ergodicity of 4> 



we have (exp tcf>)(A±) > for all large t > 0. But then (g2|) contradicts (18). So we must have 
either A + = or A_ = and by replacing ^4 by —A if necessary we may without restriction 
assume that A = \A\ = A + . Again by ergodicity we have exp tr p (cj))A = (expt(j))(A) > 
for all large t such that A > 0. The proof of non-degeneracy is now easy. Assume there 
are two linearly independent eigenvectors A\ > and A2 > 0, which we may normalize 
to 1 1^4.1 ||p = I I-A2I |p = 1. Then 7^ A\ — A 2 is also an eigenvector and by our previous 
discussion we must have either A\ > A2 or A2 > A\. This would imply that ||-Ai|| p > H^bllp 
or I \Ai ||p > I |^4i ||p respectively, again an easy consequence of @ and its proof. This is a 
contradiction, so we must have A\ = A2, thus concluding the proof of the theorem. 

4. Existence of eigenvalues 

In this section we want in particular to establish conditions for positive maps <f> in J p , which 
are sufficient to show that r p (<p) is an eigenvalue. We start with a preparation. For any A G 
B{TL) we write A = [7a|.A| with unique Ua for its polar decomposition. More precisely Ua 
is isometric on Ran|^4| = (KerA)- 1 , i.e. U* A U a'P = <p for ip € Ran|^4| = (KerA)- 1 , and 
U A <p = for tp G (Ran lA]) 1 - = KerA. 

The next two lemmas replace the Kadison-Schwarz inequality (see e.g. the second reference 
in [|TJ] and [ p4[ [| Q for other early references) used in the context of C*-algebras in [[7|, [l4j] 
and in the context of von Neumann algebras in [Q]. Since the spaces Jp, p < 00 are not algebras 
these two lemmas will provide the appropriate substitute. 

Lemma 4.1. Let the map 4> m Jp, 1 < p < 00 be 2-positive. Then for all <p,(p' £ TL and all 
A(£j p 

(23) \(<f/,<f>(A)<p)\ 2 < ^A{\A\)v)&^°^u A {\MW)- 



Also if not all terms in (23) are vanishing, then there is equality if and only if both relations 



<t>{A) tf>otfnj A (\A\) )\ -WA(A)tp) <p< 



<K|A|) <f>(A*) \f (<p, <j)(A)*ip')ip 
<f>(A) <t>o^ A (\A\) ) \ _ (^<t>(\A\)<p)<p' 



(24) ( ^JIP m t (A *L^ )( ^'A°fuj{\A\yp') y j =Q 

and 
(25) 
hold. 



For A = A* estimate (23) follows from (||), since then U\ commutes with A such that 

U A \A\U A = AU\ = U A *A = U\U A \A\ = \A\. 



Proof. Here and in the proof of the next lemma we will mimic and extend Exercise 3.4 in [26], 
page 39 in the present context. Set A x = \A\ 1 / 2 , A 2 = \A\ l l 2 U A such that A\ = \A\, A\Ay = 
A, A\A 2 = U A \A\U\ = 4>u A {\A\) are all in J p . Consider 

(26) ( Al A2 Y( Al A *]-(W A * )>0 

{ } { ) { )-{ A <p UA (\A\) 
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in TL © TL = TL ® C 2 , which defines an element in C P (TL © TL). By assumption 

Consider the following linear transformation T on C 2 (equipped with the natural scalar product) 
given as the 2x2 matrix 

(<p,<K\a\)<p) <y.0(A)V) 

(<f/,<l>(A)<p) (ip'^o4> UA (\A\)ip') 



(28) T 



such that for all a±, a 2 € 



ai \ { m \< _ ,{ ai ip \ ( 4>{\A\) (f>(A*) \ ( anp 



(29) \a 2 )> T {a 2 P- { {a 2V ' J' \<j>{A) <P o <P Ua (\A\) ) { a 2 p> P ~ ° 

with the obvious notation for the scalar product in TL © TL. A linear transformation in C 2 is > 
if and only if its trace and its determinant are both > 0. Hence the bound (|23|) follows. This 
discussion also easily gives the last claim in Lemma 4.1. Indeed, if T ^ 0, its determinant is 
equal to zero if and only if there is a non-zero eigenvector of T with eigenvalue zero and this is 
the case if and only if equality in ( |2~3| ) holds and then such an eigenvector of T is given either as 

-{<f/,<l>(A)<p) 

or as 

(¥>,#A)V> 
-(<p,<K\A\)<p) 

More precisely, if all matrix elements of T are non-vanishing then these two eigenvectors have 
non-zero entries and are proportional. If T ^ with det T = and if one diagonal element 
is vanishing, such that also the off diagonal elements are vanishing, then exactly one of these 
vectors is the null vector. In view of (|^) this concludes the proof of Lemma 4.1. □ 

Lemma 4.2. Let the map 4> in <J p , 1 < p < oo be 2-positive. Then the bound 

(30) MA)\\ p <U(\A\)\\y 2 ||0o^(|A|)||J/ a 

holds for all A £ J p . In particular A £ Ker<j) if \A\ E Ker<p or \ A\ € Ker<p o (j>u A - 
Proof. Choosing ip' = E/^m) <p in ( j23"l ) we have the estimate 

(31) {<PMA)\<P) < (v,H\A\)^) 1/2 (U^ A) vA°4u A (\A\)U^ A) v) 1/2 

for all ip £TL. Thus ( 0O| ) follows trivially from ( plj ) in case p = oo. So let 1 < p < oo. Take y>„ 
to be an orthonormal basis in TL diagonalizing |</>(^4)| with eigenvalues p n > 0. Then we obtain 

maw p = 52^ n = j2(<pn,\m\?n) p 

n n 

Y,(<Pn,H\A\)^n) P/2 (U HA) p n , 0o^(|^|)^ )¥?n )P/ 2 
n 

(32) < ||<A(|A|)||^ 2 ||0o^ A (|^|)||^/ 2 . 

Here we have used ( j3"T| ) for the first inequality. The second inequality is a consequence of 
Schwarz inequality and (^). In what follows this will become important, since we will encounter 
the situation when all inequalities are actually equalities. Also we have used the fact that with 
ip n being an orthonormal basis also Umj^Wu is a possibly incomplete set of orthonormal vectors 
and which is responsible for the last inequality. □ 



< 



< 
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A G C is called a peripheral eigenvalue of the positive map <fi in J p if A and hence also A is an 
eigenvalue with |A| = r p ((/)). 

Theorem 4.1. Let the map 4> in J p , 1 < p < oo be 2-positive with r p ((j)) = \\4>\\ p . Assume that 
A G C is a peripheral eigenvalue with eigenvector A. Then r p (<f)) is also an eigenvalue with 
eigenvector \A\. In particular, if (j> is 2-positive and compact with r p {(j)) = \ \4>\\ p , then r p (4>) is 
an eigenvalue. 

Here and in contrast to the previous discussions the case p = oo is not covered by the present 
discussion. Since J p=0 o is a C*-algebra, this case, however, is covered by the discussion in 
[36, 14]. Then and similarly in the context of von Neumann algebras one can actually prove 
more. In fact, if the spectrum is rescaled via A — > X/r p (cj)), then the peripheral eigenvalues form 
a discrete sub-group of the unit circle group, see ^ and the first reference in JI4Q. 

Proof. We want to prove that |^4| ^ is also an eigenvector with eigenvalue |A| = r p {cj)), i.e. 
^>(\A\) = r p (<j))\A\. Since \<j){A)\ = \XA\ = r p {<j))\A\ (such that also U^ A) = (X/\\\)U A and 
hence 4>u^ A) = 4>U A by tne uniqueness of U A ) it suffices to show that 

(33) 4>(\A\)ip n = H n (p n 

holds for all n, where the (p n are as in the previous lemma. Also by the previous lemma we have 

r p mA\\ p = \\\A\\ P = \\<KA)\\ P = \\\<I>(A)\\\ P 
< ||0(|A|)||i/ 2 ||^o^ A (|^|)||V2 

(34) < \\<f>\\ P \\A\\ p . 

Here we have used the fact that ||0e/ a || p < 1. By the assumption r p {4>) = \\<j>\\ p we must have 
equality. In particular this implies that 

(35) lll^)lllp = ll^l)li; /2 |bo^A(l^l)lli /2 - 

Inspection of the proof of Lemma 4.2 shows that we must have equality in (|1|) when tp = <p n 
for all n. The relations (Q) and (25 ) in Lemma 4. 1 then imply that we must have 

(36) (¥4^°^(|A|y n )0(|A|V n = ((p' n , 4>(A)ip n )4>(A)*ip f n 
and 

(37) {VnM\^\)<Pn)4> <hA\ A \W n = ^nA(A)*ip' n )c/>(A)<p n 

for all n with ip' n = U^rj^fn- Assume n to be such that p, n > 0. Since tp n € Ran |0(^4)| we 
have U^ A j ip' n = cp n and 4>(A)*ip' n = \4>(A)\ip n = fi n ip n for all such n. Similarly <j>(A)<p n = 
U ( f,(A)\4 > (A)\(pn = Hn^Pn- Hence we may rewrite ( |3^ ) and ( ^7| ) as 

(38) {ip' n ,(t>o(t>u A {\A\WM\A\)^n = l&<p n 
and 

(39) (Vn,0(|A|K)^o^(|A|V n =nlfff. 



n 



whenever \x n > 0. 

Furthermore we recall that we made use of Schwarz inequality in the estimate (32). For this 
to be an equality we claim that we actually must have 

(40) {<p n ,<l>(\A\)<p n ) = (^^o^d^K) 

for all n. Indeed, the Schwarz inequality in this context is an equality if and only if both sides of 
(40) are proportional with a common proportionality factor for all n. By the equality (35) this 
proportionality factor has to be equal to 1, thus proving (pQ|). 

If p n = then ip' n = and hence by ( |40| ) c/)(A)<p n = 0. Thus we have established (33) in the 
case when p n = 0. If p n > then the left hand sides of (38) and (39) are both non- vanishing 



PERRON-FROBENIUS THEORY FOR POSITIVE MAPS ON TRACE IDEALS 



11 



since then ip' n / 0. Therefore ip n is an eigenvector of i^(|^4|) with eigenvalue 

(<p n ,<l>(\A\)<p n ) - 



(cp' n ,<po(f) UA (\A\)ip' n ) 
2 



(<p n ,(j)(\A\)(p n ) 

This gives ([53|) in the case when fi n > 0, thus concluding the proof of the theorem, since for 
compact operators the spectrum is pure point. □ 

The use of (|7|) in the second inequality in (32) does not give any additional mileage when p / 1 
since the cp n are eigenvectors of <^(|^4|) and the tp' n eigenvectors of 4> o <ftu A and hence (fy gives 
an equality for the case at hand. 

The strategy of the proof was motivated by the following observation in the classical context. 
With the notation introduced after (||) assume that Sz = Xz. This gives |A| \z\ = \Xz_\ = \Sz\ < 
S\z\. Taking its usual Hilbert space norm in C ra we obtain |A| 1 1 \z\ \\ < \\S\z\ \\ < \\S\\\\\z\\\. 
So if |A| = 1 15|| we must have S\z\ = \\S\\ \z\. 

The next result extends Theorem 3.1. 

Theorem 4.2. Let the map <fi in J p , 1 < p < oo be 2-positive and ergodic with r p {(p) = 
\\4>\\ p - Assume that X € C is a peripheral eigenvalue. Then X is a simple eigenvalue and the 
corresponding eigenvector A is a normal operator. 

Proof. Assume there are two eigenvectors A\ and A2 in J p with eigenvalue A, normalized to 
||j4i|| p = || I Ai| ||p = I \A2 ||p = || \A% 1 ||p = 1. By Theorem 3.1, its proof and by Theorem 4.1 
we have that < = IA2I is an eigenvector with eigenvalue r p (</>). In particular both Ua 1 
and Ua 2 w& unitaries. Since A\ — exp it A2 is also an eigenvector for the same eigenvalue for 
all < r < 2ir by the same argument we must have 

(41) I — expir A2 1 = \ \A\ — exp it A2\\ p \Ai\. 

The function r h-> \\Ai — expir A2\\ p is easily seen to be continuous. Hence there is to such 
that 

(42) ||^i — expiro A^Wp = min \\A\ — expir A2IL. 

0<T<2lT 



Take the square of ( |41| ) and use ( |42| ) to obtain 

|^4i — expiro A2 1 2 < \A± — expir A2I 



Due to the relation Ran | A\ \ = H this gives 

\\(Ua 1 - expiro U A2 )(p\\ 2 < IK^Ai - expfr Ua 2 )v\\ 2 
which written out in turn gives 

(43) K (e lT (U Al <p, U A2 if}) < K (f*(U Al tp, Ua 2 v)) 

for all < r < 2ir and all ip G Tt. We distinguish two possible cases. First if A% — exp iro A2 = 
and hence Ua x — exp iro U a 2 = there is nothing to prove. Otherwise U = exp iro Ua 2 
is a unitary operator 7^ I and hence by the spectral theorem for unitary operators there is ip G ri 
such that 3ft( (if, U ip)) < \(ip,U<p}\. Choose r such that exp i(r — ro )(</?, Uip) = \{tp, U(p}\. 
This contradicts (p3|), thus concluding the proof of the first part of the theorem. As for the last 
part XJ a is unitary, as already noted. Relation ( |39| ) shows that </>(/ a (|t4|) is also an eigenvector 
of (f) with eigenvalue r p {4>). Hence ^^(l^l) = \A\ by uniqueness and normalization. So Ua 
commutes with \A\ and therefore A is normal with U\ = Ua* and \A*\ = \A\, where we recall 
that A* is an eigenvector with eigenvalue A. □ 

The next result is of relevance in quantum physics. We recall that for completely positive, 
trace preserving maps in J\ one has r\{(j)) = ||0||i = 1. So we immediately obtain the 
following existence result for Perron-Frobenius vectors. 

Corollary 4.1. Let the map 4> in J\ be completely positive, trace preserving, compact and er- 
godic. Then there is a unique density matrix p > which is left invariant under (f). 



12 



R. SCHRADER 



A converse of this result is given by Theorem 5.1 below. 

The following result is of the type needed in the standard context of Monte Carlo simulations 
in lattice models of statistical mechanics. It states that, starting from any non-negative initial 
configuration, by iterations of an ergodic up-dating procedure (heat bath method) one eventually 
arrives at the desired equilibrium configuration which typically is given by a Gibbs distribution 
and which by construction of the up-dating is a Perron-Frobenius vector for the up-dating pro- 
cedure. 

Corollary 4.2. Let the map <p in J 2 be 2-positive, self adjoint, compact and ergodic with Perron- 
Frobenius vector A > normalized to \\A\\2 = 1. Then for any B € J 2 

(44) s - lim exp(tf>- M\ 2 ))(B) = (A,B) 2 A. 

t^oo 

If < B ^ then the right hand side of @ is non-vanishing. It would be interesting to find 
an analogous formulation when p 7^ 2. 

Proof. We use familiar arguments. Let the Aj form an orthonormal basis of eigenvectors for (j) 
arranged such that A\ = A and <x/+i < aj. Here the afs are the eigenvalues, i.e. (j>(Aj) = 
CFjAj. Since a\ = \\(t>\\ 2 is a simple eigenvalue, o\ > a 2 . By Plancherel's theorem we have 
B = J2j bjAj with bj = {A h B) 2 and J2j = \\B\\%. Obviously 

exp(^- \\<t>\\ 2 ))(B) = Y J b 3 e t ^~^ ) A, 

j 

holds such that 

(45) ||(exp(t(^- ||^|| 2 )))(B) - (A,B) 2 A\\ 2 < e^-^WB^ 

and the claim follows. □ 



By a slight modification of the proof using the spectral representation of one may replace 
the compactness condition by the condition that ||0||2 is an isolated eigenvalue, i.e. \\4>\\2 is 
separated by a distance m 2 (the "mass gap" in physical language) from the remainder of the 
spectrum of (j), which otherwise may be arbitrary. Then the estimate (45) is still valid with 
<7i — o"2 being replaced by m 2 . 



5. Examples 

In this section we will provide non-trivial examples when dimH = 00. The issue is to show 
that all 4 conditions in Corollary 4. 1 can be fulfilled simultaneously. In the finite dimensional 
case the compactness criterion is automatically satisfied and then it is easy to construct non- 
trivial examples. In general the existence of completely positive, trace preserving, compact 
maps in J\ is clarified by 

Lemma 5.1. Let a = {ai}i^ be such that YlieN a i a i = I an d all ai € B(7i) are compact. 
Then the map 4> a in J\ is compact. 

Proof. We recall that cj) a is the norm limit of </> Qjv as ./V — > 00. Since the space of compact 
operators is closed w.r.t. the norm topology it suffices to prove that all Qjv are compact. By 
assumption to each i there is a sequence of finite rank operators in BTL) such that 

lim \\a.i — otikW = li m \\ol* — oi£k\\ = 0- 

fe— >oo ' k— >oo 

Since a* k is of finite rank <fia Nk with ajvfc = ( a l,k> a N,k) is also of finite rank and in 
particular compact. Using the a priori estimate | |ABC| |i < \\A\ \ \\B\\i ||C|| and a standard e/3 
argument gives 

Urn H^-^JIi =0. 

□ 
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We are now prepared to provide non-trivial examples, which resulted from a discussion with D. 
Buchholz. Let ipi, i € N be any orthonormal basis in TL and c- lk > any set of numbers which 
satisfy the condition J2i c< i k = ^ f° r a ^ ^. m tne Dirac notation we define the following family 

a = {aik}i<i y k<oo of rank 1 operators 

«ifc = c ik \^i)(ip k \, 

such that 

o4 = <Hk\il>k)(il>i\, 

In other words is the map (p t— » Ci k (tp k , We set 

0«(-4) = ^2 aikAa* k 

i.k 

By assumption we have 

J2a* k a ik = J2cf k iV'fcXV'fcl = \^k)(ipk\ = I, 

so ^>q, is completely positive, compact and trace preserving. We claim that it is also positivity 
improving by Lemma 3.3. Indeed, given ^ <p € TL, assume there is ^ ifo G TL such that 

= ((p ,a* k ip) = c ik ((p ,ip k )(ipi,ip) 

holds for all i and k. Since by assumption Ci k > this contradicts the fact that the ipi form an 
orthonormal basis. So this (f> a satisfies all conditions of Corollary 4.1. The resulting Perron- 
Frobenius vector is given as the density matrix 

(46) p = ^2pi\ipi){ipi\, Pi > 0,^2 Pi = 1, 

i i 

where the pi form the components of the Perron-Frobenius vector with eigenvalue 1 for the 
matrix cf k , i.e. J2k c lkPk = Pi- As a special case consider the choice eg. = q > for all k with 
Ejcf = 1. Then p { = cf. 

This discussion also gives the following converse to Corollary 4.1, which in its standard 
formulation forms one of the starting points for Monte Carlo simulations in lattice models of 
statistical mechanics and (euclidean) quantum field theory including lattice gauge theories. 

Theorem 5.1. Given an arbitrary density matrix p > 0, there is a completely positive, trace 
preserving, compact and ergodic map <j) in J\ for which p is the resulting Perron-Frobenius 
vector. In addition <fr may be chosen to be idempotent (i.e. <fi satisfies the characteristic equation 
cj) o (1 - <f>) = 0) and o-i ((f)) = {0, 1}. 

Proof. Take the ^'s to form an orthonormal basis in which p is diagonal, i.e. such that p 



1/2 

takes the form fl46| ) and choose cf) = <pg_ with a as above with Ci k = p i . Then we have 
4>(A) = Tr A ■ p for all A € J\, such that is idempotent. Also the kernel of consists of all 
elements which may be written in the form A — Tr A ■ p. Moreover let A ^ 0, 1. Then for given 
A' the equation (A — 4>){A) = A 1 has a solution in A of the form 



□ 



As in the corresponding ( finite dimensional) standard context cj) with prescribed Perron- 
Frobenius vector is of course not unique. We provide another example which in spirit is 
much closer to the usual "local" up-grading procedure in Monte Carlo simulations. Set -p v ^ = 
p u+i /(p i+ i +pi + pi-i) for v = 0, ±1 and with the convention p = 0. Then J2u=o,±i Vv,% = 1 
for all i and p Ui i > for all v and i unless u = — 1 and i = 1. Set 

T v,i =pl,i \i>u+i){lpi\ 

and define 4>t as 

v,i 
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An easy calculation shows that cfyp is completely positive, trace preserving, compact and ergodic. 

Acknowledgments The author has profited from discussions with D. Buchholz, B. Kummerer, 
M.B. Ruskai and E. St0rmer. 
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